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OPERATING  CHARACTERISTICS  FOR  CLIPPING  OF  IN-PHASE 
AND  QUADRATURE  COMPONENTS  OF  INPUT  AND/OR 
REFERENCE  OF  NARROWBAND  CORRELATOR 

INTRODUCTION 

To  minimize  the  amount  of  data  processing,  telemetry  bandwidth,  execution 
time,  and  storage,  hard  clipping  is  frequently  employed  in  signal  processing 
hardware.  In  particular,  for  cross-correlation  of  a  received  input  waveform 
with  a  local  reference,  several  alternatives  exist  for  inclusion  of  clipping. 
Either  the  input  or  the  reference  or  both  could  be  clipped.  For  a  narrowband 
system,  where  the  input  is  complex  demodulated,  the  further  option  of 
employing  clipping  in  both  the  in-phase  and  quadrature  channels,  of  the  input 
as  well  as  reference,  is  available. 

Here  we  will  investigate  all  the  various  possibilities  of  including 
clipping  at  the  input  and/or  reference  level,  for  both  baseband  (real)  and 
narrowband  (complex)  signal  processing  systems.  In  this  manner,  we  will 
ascertain  whether  inclusion  of  the  additional  clippers  in  the  quadrature 
channels  ameliorates  the  degradation  associated  with  the  usual  case  of  a 
single  clipper  operating  on  one  real  input  process. 
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For  completeness,  we  will  present  the  analysis  and  simulation  for  both 
the  baseband  correlator  as  well  as  the  narrowband  correlator.  This  will  serve 
as  verification  of  the  analysis  technique  and  afford  a  ready  comparison  of 
both  types  of  correlators. 

No  definitions  of  output  signal-to-noise  ratio  criteria  are  employed 
here.  Rather,  we  evaluate  the  system  output  detection  probability  and 
false  alarm  probability  Pp  in  terms  of  the  input  signal-to-noise  ratio, 
which  is  a  well  defined  quantity.  We  can  then  make  a  direct  comparison 
between  systems,  of  the  required  input  signal-to-noise  ratios,  in  order  to 
achieve  some  specified  common  performance  level  Pr,,Pr- 
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PROBABILITIES' FOR  LINEAR-INPUT  CHANNEL 


In  this  subsection,  we  presume  that  the  nonlinearity  g  in  the  input 


channel  of  figure  1  is  absent;  that  is,  from  (2)  and  (3),  The 


output  of  the  correlator  is  then,  for  signal  present. 


'  ■  2  “k'i=  ■  ^  "k^k  '  2  “k'^k  "  "k’  ' 


(5) 


k  k  k 

where  ^  denotes  the  sum  from  k=l  to  K.  It  is  important  to  observe  that 

K 

signal  waveform  [sj^]  and  reference  [wj^j  are  completely  general  at  this 
point.  That  is,  we  have  not  restricted  consideration  to  the  examples  of  (4) 
yet . 


Since  process  [n  ]  is  Gaussian,  the  random  variable  z  in  (5)  is  also 


Gaussian.  It  has  mean  and  variance 


=  2  ^ 


k^k 


=  a 


n 


(6) 


(The  latter  relation  also  holds  true  when  the  signal  is  absent.)  The 
detection  probability  is  then,  for  threshold  T, 

,2 


oC 


Pp  ^  Prob(z  >  T)  - 


du 


17 

j  (2ir) 


exp 


(u  -  m^) 


2a 


(7) 


where  normal  probability  integral 


dt  {2n)  exp  -t  /2  ; 


dt  0( t ) 


:8) 


-  00 


f' d r,  •'d  '^d 'T. ^wMWdwvwvw-iwuf  x* '•(.■V.* HT •*."  V H."  V .^TTrjVLV^ V"3ri,V.^ '.^^.^.'V'ji 
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and  deflection  parameter 


The  false  alarm  probability  is  obtained  by  setting  the  signal  |S|^^ 
equal  to  zero  in  (9)  and  (7): 

Pp  =  $(-T/a^)  .  (10) 

Combining  (7)  and  (10),  we  have  the  exact  performance  (receiver  operating) 
characteristics  in  the  compact  form 

Pp  =  $(d  ^l(Pp))  .  (11) 

where  $  is  the  inverse  function  to  ^  of  (8).  Thus,  the  single  parameter  d  in 
(9)  and  (11)  is  a  complete  descriptor  of  performance  for  a  linear  input 
channel  in  figure  1.  Furthermore,  [Sj^j  and  [wjJ  are  completely  general  in 
(9).  Observe  that  the  absolute  scale  of  the  reference  cancels  out  in  (9). 

Matched  Reference 


If  the  local  reference  is  matched  to  the  received  signal,  then 
w,  =  As,  as  in  (4),  and  (9)  yields  the  matched  deflection  parameter  value 


where  E  is  the  received  signal  energy.  Furthermore,  this  is  the  maximum 
possible  value  of  d  in  (9)  by  choice  of  local  reference  Thus,  the 

optimum  performance  of  the  linear-input  correlator  depends  on  the  total 
received  signal  energy,  and  not  on  its  f ract iona 1 i zat ion  into  individual 
samples 

Clipped  Reference 

For  this  case,  (4)  yields  =  A  sgnfSj^),  and  (9)  specializes  to  value 


n 


This  quantity  depends  on  the  specific  fractional i zation  of  the  received  signal 
energy  into  components  {s.]  and  is,  therefore,  example-dependent.  Its 
maximum  value  is  again  E  as  in  (12),  but  only  if  all  components  are 

equal  in  magnitude.  We  will  consider  two  particular  signal  examples  in  the 
rest  of  this  section  on  the  baseband  correlator. 

EXAMPLES 

Example  1,  Cosine  Wave: 

S|^  =  a  cos  ^  2Trk/K  for  1  £  k  <  K  ;  a  >  0  .  (14) 


Example  2,  Constant; 
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For  these  two  examples,  the  values  of  the  matched  deflection  parameter  in 
(12)  become,  respectively. 


m 


(2)  ^  ^1/2  a_ 


m 


(16) 


On  the  other  hand,  for  a  clipped  reference,  (13)  yields  corresponding  values 


.(1)  _L 


r  ^1/2 


^  |cos(2irk/K)| 


7.201  a/a^  for  K  =  128) 


k=l 


n  1 .707  a/a  for  K  =  8 
n 


d(2)  ^  ^1/2  ^ 


for  a  1 1  K  . 


(17) 


n 


Since  d  is  a  complete  descriptor  of  performance  for  the  linear-input 
correlator,  as  given  in  (11),  it  is  seen  that  for  the  cosine  signal  example  1, 
the  clipped  reference  requires 


-20  log 


20  log(7. 201/8)  =  .91  dB 

20  log(l .707/2)  =  1  .38  dB 


for  K  =  128 
for  K  =  8 


(18) 


additional  input  signal -to-noise  ratio  relative  to  the  matched  reference.  On 
the  other  hand,  constant  signal  example  2  requires  a  0  dB  difference,  as  seen 
by  reference  to  (16)  and  (17).  In  general,  (12)  and  (13)  reveal  a  dB 
difference,  due  to  clipping  the  reference,  of 


10  log  K 


(19) 
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PROBABILITIES  FOR  NONLINEAR  INPUT  CHANNEL 


Here,  the  nonlinearity  g  in  the  input  channel  of  figure  1  is  present; 
that  is,  V|^  =  g(X|^)  as  in  (3).  Not  yet  specializing  to  the  clipper  of 
(2),  we  have  system  output 


^  =  ,2  “k  3^*k^  =  ^  ^  • 


(20) 


where  signal  [sjJ  and  reference  [wj^j  are  also  general. 


The  exact  evaluation  of  the  distribution  of  random  variable  z  in  (20)  is 
difficult;  accordingly  we  limit  consideration  to  the  case  where  K  »  1, 
meaning  that  z  is  approximately  Gaussian.  The  p-th  moment  of  the  output  of 
nonlinearity  g,  conditioned  on  input  signal  value  S|^,  is 


*  "k* 


(2*) 


1/2 


a 

n 


g^(s^,  n)  = 


dx  (S(x)  g  (s. 


(21) 


where  we  used  the  Gaussian  character  of  noise  {n|^]  and  (B).  Accordingly, 
the  mean  and  variance  of  random  variable  z  in  (20)  are 


^  ^  '^k  ^l^'k)  ' 

k 


•>1  ■ 


'  "k>) 


^2^^k^ 


r2. 


n. 


(22) 


using  the  independent  identically  distributed  property  of  noise  fn  ]  . 
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Then  in  a  manner  similar  to  (7),  the  detection  and  false  alarm 
probabilities  are  given  approximately  by 


D 


“f-T 


I'm  -  T’ 
zo 

®zo 


(23) 


where  T  is  the  threshold,  and  the  sub  o  designates  setting  the  signal  (s^ 
to  zero  in  (21)  and  (22).  Upon  elimination  of  threshold  T,  the  receiver 
operating  characteristic  is  governed  by  the  relation 


P 


0 


(24) 


where 


is  the  appropriate  deflection  parameter  in  this  case  of  a  nonlinearly 
distorted  input  channel.  However,  a  is  not  a  complete  descriptor  of 
performance,  since  output  z  in  (20)  is  not  precisely  Gaussian;  however,  for 
large  numbers  of  samples,  K,  the  receiver  operating  characteristic  furnished 
by  (24)  should  be  a  good  approximation. 


PROBABILITIES  FOR  CLI PPEO-INPUT  CORRELATOR 

We  now  specialize  the  nonlinearity  in  figure  1  and  (20)  to  the  clipper  of 

(2).  Then  (21)  yields  the  first  two  moments 
10 


TR  7757 


G^(s)  =  Jdx  0(x)  sgn(s  +  o^x)  =  ^  • 

62(5)  =  Jdx  0(x)  sgn^(s  +  o^x)  =  1  . 


(26) 


Substitution  of  these  results  in  (22)  then  yields  mean  and  variance 


Upon  setting  the  signal  {s|j}  to  zero  in  (27).  there  follows 


m 


zo 


=  0 


(27) 


(28) 


Utilization  of  (27)  and  (28)  in  (24)  and  (25)  yields  the  results  for  the 
clipped-input  correlator,  namely. 
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Approxifnation  for  Small  Signal -to-Noise  Ratio 


If  the  input  signal-to-noise  ratio  in  figure  1  is  small,  that  is, 


«  1  for  all  k 


(31) 


there  follows  from  (8),  the  approximation 


$ 


=  \  ^  (2ir) 


-1/2 


Then  (29)  simplifies  to 


(32) 


(33) 


the  last  relation  following  directly  from  (9).  And  since  (30)  now  reduces  to 

o  /o  21,  (24)  shows  that  A  now  is  the  sole  parameter  describing  the 
zo  z 

receiver  operating  characteristic,  under  the  assumptions  of  large  K,  number  of 
samples,  and  small  input  signal-to-noise  ratio. 


Equation  (33)  reveals  that  the  clipped-input  system  of  figure  1  is 
degraded  relative  to  the  linear  system  by 

ri  / 

20  log  ^  =  20  log(l)  =  1.96  dB  .  (34) 

regardless  of  the  particular  signal  {sj  and  regardless  of  the  reference 
employed.  Thus,  clipping  the  input  to  the  correlator  in  figure  1 
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causes  a  loss  of  1.96  dB  when  K  is  large  and  the  input  signal -to-noi se  ratio 
is  small,  no  matter  what  signal  and  reference  are  used. 

Matched  Reference 


The  deflection  parameter  a  In  (33)  Is  maximized  by  choosing  a  matched 
reference,  namely,  w^^  =  As^^,  thereby  yielding 


The  last  identity  follows  from  (12).  Again,  the  performance  of  the 
c 1 i pped-i nput  correlator  depends  only  on  the  total  received  signal  energy,  and 
not  on  its  fractional i zation  into  components  {S|^]  • 

Clipped  Reference 

For  the  clipped  reference,  with  Wj^  =  A  sgn(Sj^)  from  (4),  (33)  reduces 
to 

“r  ■  <3'') 

'’n  k 

the  last  relation  following  from  (13).  The  loss  relative  to  the  linear-input 
c  I  ipped -ref erence  correlator  is  again  1.96  dB,  since  (35)  and  (36)  are  special 
cases  of  (33)  and  ( 34 ) . 


13 
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EXAMPLE 


For  the  cosine  signal  already  considered  in  (14),  the  deflection 
parameters  in  (35)  and  (36)  become 

/IX  /•x\l/2/,/\1/2  ,  /o\l/2 


x(1)  = 


5.746  a/o 


1.362  a/o. 


for  K  =  128 


for  K  =  8 


see  (16)  and  (17).  The  first  line  reveals  a  1.96  dB  loss  relative  to  the 
linear-input  correlator,  both  with  matched  references.  The  second  line 
indicates  a  loss  -dependent  on  the  particular  number  of  samples,  K. 


GRAPHICAL  RESULTS 


In  figure  2,  the  receiver  operating  characteristics  (ROC)  for  the 
linear-input  correlator,  as  given  by  (11),  are  drawn  in  dotted  lines,  for  the 
range  (.001,  .8)  in  false  alarm  probability  and  (.001,  .999)  in  detection 

probability.  Superposed  as  the  jagged  solid  lines  are  the  results  of  two 

* 

simulations,  each  employing  30,000  trials  ,  for  the  linear-input  correlator 

with  K  =  128  samples  and  for  the  cosine  signal  example  of  (14).  For  case  A, 

a/o  =  .5,  whereas  for  case  8,  a/o  ^  .25.  Reference  to  (16)  reveals  that 
n  n 

these  values  correspond  to  deflections  d^^ ^  =  4  and  d|^^ ^  =  2,  respectively. 


*A  sample  program  is  listed  in  the  appendix. 


Prcbability  of  False  Rlarrs 

2.  ROC  for  Baseband,  Linear,  K=128 


and  in  fact,  the  simulations  overlay  these  d-values  very  closely  in  figure  2, 
over  the  entire  range  shown.  This  agreement  establishes  the  degree  of 
confidence  to  be  expected  in  the  simulation  results  to  follow. 


In  figure  3,  the  value  of  K  is  decreased  to  8,  while  the  values  of  a/o 

n 

are  increased  by  a  factor  of  4.  According  to  (16),  this  results  in  precisely 
the  same  d-values  as  above,  and  again,  the  simulations  overlap  the 
corresponding  analytic  results,  even  for  this  small  value  of  K,  as  they 
should,  since  the  correlator  is  linear. 


In  figure  4,  K  is  increased  to  128,  but  the  reference  is  now  clipped. 
According  to  the  results  in  (17),  we  have 


(3.60 

for 

a/<,„  = 

.5 

|l  .80 

for 

a/a„  • 

.25 

for  K  =  128  . 


These  values  are  borne  out  by  the  simulation  for  the  cosine  signal  example  in 
figure  4. 


The  corresponding  simulation  results  in  figure  5  overlay  these  d  values  over 
the  entire  range  plotted. 
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Figure  4.  ROC  for  Baseband,  Clip  Reference,  K-128 
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When  the  input  is  dipped,  instead  of  the  reference  which  is  matched 
according  to  (35),  the  pertinent  equation  is  the  upper  line  of  (37).  There 
follows 


r3.i9 

for  a/a  = 
n 

.5 

|l.60 

for  a/a^  = 

.25 

for  K  =  128  . 


Figure  6  correlates  these  values  very  well. 


If  K  is  decreased  to  8,  while  a/o  is  increased  to  2  and  1, 

n 


respectively,  the  same  values  result  from  the  use  of  (37),  as  given  in  (40). 
The  simulation  results  in  figure  7  reveal  that  the  actual  performance  does  not 
meet  those  deflection  values  3.19  and  1.60,  except  at  the  upper  right  ends  of 
the  curves.  Two  obvious  reasons  for  the  discrepancy  are  that  K  =  8  is  not 
large  enough  to  rely  on  the  Gaussian  approximation,  and  the  input 


signal-to-noise  ratios,  a/o^  =  2  or  1 ,  are  not  small,  as  assumed  in  (31) 


et  seq. 


In  an  effort  to  circumvent  the  small  input  signal-to-noise  ratio 
assumption  made  in  (31),  we  returned  to  the  more  general  results  (29)  and  (30) 
and  employed  them  in  receiver  operating  characteristic  (24).  The  values 


calculated  for  the  matched  reference  w.  =  As,  were 

k  k 


4.932  for  a/a  =  2 
n 


1.785  for  a/o^  =  1 


for  K  ^  8  , 


(2.298  for  a/o 
<3  n 

zo  J 


1.258  for  a/o.  ^  1 


for  K  -  8 


Probability  of  Detection 
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These  results  are  superposed  in  figure  7  with  the  label  (24,  29,  30),  and  are 
seen  to  be  poorer  than  those  obtained  via  the  low  input  signal -to-noise  ratio 
results.  However,  this  is  felt  to  be  a  fortuitous  circumstance. 

When  both  the  input  and  the  reference  are  clipped,  the  simulation  results 
for  K  =  128  are  displayed  in  figure  8.  The  corresponding  theory  is  furnished 
by  (36)  and  yields  the  results  for  ^  already  listed  in  (37).  They 
become 


for  K  =  128  ,  (42) 


p.873 

for 

a/o  = 
n 

.5 

|l  .436 

for 

a/On  = 

.25 

in  excellent  agreement  with  figure  8. 


For  the  alternative  situation  with  K  =  8  and  a/o  =  2  or  1 ,  the 

n 


appropriate  values  follow  from  (37)  as 

r2.724  f 

.(1)  ,  } 


2.12^  for  a/o  =  2 
n 


1.362  for  a/o^  =  1 


for  K  =  8  . 


These  values  do  not  furnish  a  good  approximation  to  the  corresponding 
simulation  results  depicted  in  figure  9,  except  at  the  upper  right  end. 
Hov;ever,  when  we  resort  to  the  more  accurate  approach  of  (24),  (29),  (30),  as 
discussed  above,  agreement  is  considerably  better,  although  the  slopes  of  the 
theoretical  curves  are  steeper  than  those  of  the  actual  simulation  results.  A 
similar  situation  occurred  in  figure  7,  although  considerably  worse  there. 


■f^r 
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NARROWBAND  CORRELATOR 


In  this  section,  the  correlator  input  and  local  reference  are  complex 
sampled  processes,  corresponding  to  the  in-phase  and  quadrature  components  of 
the  narrowband  processes;  see  figure  10.  The  double  arrows  denote  a  pair  of 
samples,  as  for  the  complex  envelope  of  a  narrowband  process.  It  is  presumed 
that  any  time  delays  or  frequency  shifts  of  the  input  signal  have  been 
compensated  for,  and  we  concentrate  on  the  effects  of  clipping,  at  various 
locations,  on  the  receiver  operating  characteristics  of  the  correlator. 


Input  Nonlinearity  Reference 


Accumulator  Magnitude 


Figure  10.  Narrowband  Correlator 


SYSTEM  DESCRIPTION 


The  input  to  the  correlator  in  figure  10  is 


n,  for  Hr, 
k  0 


exp(  ie)  n|^  f  or  H, 


where  [sj  is  a  deterministic  known  complex  waveform,  and  e  is  a  random 
variable.  The  phase  shift  e  is  independent  of  k;  that  is,  the  phase  shift  is 
constant  but  unknown  over  the  observation  interval  K. 
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Additive  noise  {n  ]  is  complex  (circularly-symmetric)  Gaussian  with 


independent  identically  distributed  components: 


"k  =  ° 


for  a  1 1  k  , 


n.  n  =0 
k  m 


*  2 

n,  n  =  2a  5. 
km  n  km 


for  all  k,  m  . 


2  2 

In  particular,  this  means  that  |n|^  =  2a^. 


Expressing  the  noise  in  terms  of  its  real  and  imaginary  parts, 


n,  =  n,  +  1  n,  .  , 
k  kr  ki 


(45)  yields  the  properties 


"kr  ■  "ki  =  “ 


n,  n,  .  =  0 
kr  ki 


for  al 1  k  , 


2  2  2 

n,  =  n,  .  =  o„ 
kr  k  1  n 


n,  n  =  n,  .  n  .  ^  n.  n  .  =  0 
kr  mr  ki  mi  kr  mi 


for  all  'k  m 


Thus,  is  the  variance  of  each  component  of  the  received  noise  fn.1 . 
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The  nonlinearity  g  in  figure  10  is  as  given  earlier  in  (2).  However, 
since  there  are  both  in-phase  and  quadrature  channels  input  to  the 
nonlinearity,  namely,  and  x^^.,  its  output  is  also  complex: 


%  ■  ’'-kr’  "  ’  ’‘"kl*  '  “kr  "  '  Vi  • 


Thus,  two  nonlinearities  are  employed,  one  for  each  channel.  When  g  is  the 
clipper  given  by  (2),  output  Vj^  takes  on  one  of  the  four  values  ±1  f  i,  a 
2-bit  representation. 


Local  reference  {w  ^  is  also  complex  deterministic  and  directly  related 


to  the  signal.  Two  examples  are 


for  linear  reference 


A  sgn(S|^^)  i  sgn(S|^^)  for  clipped  reference 


Again,  two  clippers  are  required  in  the  reference  channel,  one  for  each 
component.  This  also  results  in  a  2-bit  representation  for  the  reference. 
Complex  scale  factor  A  is  irrelevant  to  the  performance  of  the  narrowband 
correlator. 


The  multiplier  output  in  figure  10  is  given  by  W|^V|^;  that  is,  a 
conjugate  is  applied  to  the  local  reference.  This  complex  quantity  is  summed 
over  the  observation  interval  of  length  K,  yielding  complex  output 


w*v.  = 


Z  f  1  z . 
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Finally,  the  magnitude  of  this  quantity. 


is  compared  with  a  threshold  for  decision  about  signal  absence  (H^)  or 
presence  (H^. 

PROBABILITIES  FOR  LINEAR-INPUT  CHANNEL 

In  this  subsection,  we  presume  that  the  nonlinearity  g  in  the  input 
channel  of  figure  10  is  absent:  that  is,  from  (2)  and  (48),  Then 

accumulator  output  z  in  (50)  becomes,  for  signal  present, 

'  =  5  Vk  =  ^  Vk  =  ^  '"k('k  +  n^)  =  S  f  N  ,  (5: 

k  k  k 

where 


S  =  exp(ie)  ^  W|^S|^  ,  N  =  .  (5 
k  k 

It  should  be  observed  that  complex  signal  waveform  [sj  and  local  reference 
[wj  are  completely  general  at  this  point.  That  is,  we  have  n^  restricted 
consideration  to  the  examples  given  in  (49)  yet. 


r,  W,'  w  w.-l 
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For  a  given  fixed  e,  random  variable  z  in  (52)  is  complex  Gaussian,  since 
the  operation  of  summation  is  linear  in  the  noise  variables  {n|^}  .  We  will 
find  the  probability  of  detection  conditioned  on  a  fixed  value  of  e.  The 
complex  random  variable  N  in  (53)  is  Gaussian  with  moments 


^ —  ir  ★  - 

>  w,  w  n,  n  =  0  , 
^  k  m  k  m 


.  5  w*w  n,n*  .  5  Iv,  2.1  . 

^  kmkm  n^:llkl  N 


Here  we  used  the  noise  properties  listed  in  (45).  The  statistics  given  in 
(54)  hold  true  whether  signal  is  present  or  not  at  the  correlator  input. 


If  we  let  N  be  expressed  in  terms  of  its  real  and  imaginary  parts, 


N  =  N  +  i  N.  , 
r  1 


the  results  of  (54)  translate  into 


N  =  N.  =  0  , 
r  1 


N^N.  =  0  ,  nJ  =  .  (56) 

Therefore,  we  can  write  the  joint  probability  density  function  of  the  real  and 
imaginary  parts  of  z  in  (52)  as 
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:,.z,)  .  (2,0^)  exp  - 


L- 


for  signal  present. 


The  detection  probability,  for  threshold  T,  is  given  by 


Pp  =  Prob(Y  >  T)  =  Prob 


>  T  =  dz^  dz.  p(z^.z.) 


if  "'r  0^1 


exp  - 


4  *  y  -  "  1^1 


=  (2ir) 


i  dp  p  j  exp|-  f-  +  ^ 

T/a  -ir 

n 


'  \  Isl ^ 

coses  +  sinaj  f  -^-y 


dp  p  exp  [  l(p^  +  d^)  I^(dp) 


Q(d,  T/P^)  , 


where  we  used  in  order;  (51),  (57),  the  substitution 


z^  =  P  COS0  ,  z.  =  p  sina  , 


[I;  8.431  3],  and  Marcum's  Q-function  [2;  (I)  for  M  ^  I],  The  parameter 


introduced  in  (58),  namely. 
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is  the  analogue  here,  in  the  narrowband  correlator,  to  the  quantity  (9) 
defined  for  the  baseband  correlator;  we  made  use  of  (53)  and  (54)  in  (60). 

Although  the  derivation  in  (54)-(58)  was  conditioned  on  a  specific  fixed 
value  of  e,  the  end  result  in  (58)  and  (60)  depends  only  on  js)  which  is 
independent  of  e.  Therefore,  the  detection  probability  of  the  narrowband 
correlator  is  independent  of  whatever  e  is,  and  (58)  is  also  the  unconditional 
detection  probability,  there  being  no  need  to  average  over  e;  its  exact 
probability  density  function  is  irrelevant. 


The  false  alarm  probability  is  obtained  from  (58)  and  (60)  by  setting  the 


signal  [s^^  to  zero,  and  using  [2,  (2)]: 


=  0(0,  T/oJ  =  exp  - 


Eliminating  T/o^  from  (58)  and  (61),  we  obtain  the  operating  characteristic 
in  the  form 


(-2JJn  Pp) 


(62) 


Thus,  deflection  parameter  d  in  (60)  is  a  complete  descriptor  of  performance 
for  the  linear-input  narrowband  correlator  in  figure  10.  Furthermore,  signal 
[S|^^  and  local  reference  [W|^]  are  completely  general  in  (60). 
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Matched  Reference 


If  the  local  reference  is  matched  to  the  received  signal,  then 
=  AS|^  as  in  (49),  and  (60)  yields  for  the  matched  deflection  parameter 
the  value 


where  E  is  the  energy  of  the  received  real  signal.  Furthermore,  this  is  the 
maximum  possible  value  of  d  in  (60)  by  choice  of  local  reference  . 

Thus,  the  optimum  performance  of  the  linear-input  narrowband  correlator 
depends  on  the  total  received  signal  energy,  and  not  on  its  fractional i zation 
into  individual  components 

Clipped  Reference 


Since  (60)  and  (62)  hold  for  any  reference,  we  can  specialize  it  to  the 
clipped  case  given  in  (49),  namely. 


"k  -  “farKSkr)  »  '  ■ 


(64) 


Substitution  in  (60)  yields 


d,  =  (2K) 


J  [sgn(s^^)  -  i  sgn(s^.)][s^^  .  i  s^.j 


(65) 


This  quantity  depends  on  the  specific  f rac t ional i zat i on  of  the  received  signal 
energy  into  components  and  is,  therefore,  example-dependent.  Its  maximum 
value  is  again  (2E)  as  in  (63),  but  only  if  all  components  [S|^] 

are  the  same  complex  constant.  The  degradation  of  d^  relative  to  d^ 
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depends  on  the  specific  signal  .  We  will  consider  two  particular  signal 
examples  in  the  rest  of  this  section  on  the  narrowband  correlator. 

EXAMPLES 


Example  1,  Phase-modulated  tone; 

s^  =  a  exp(ii8,^).  =  2irk/K  for  1  <  k  <  K;  a  >  0  .  (66) 

Example  2,  Pure  tone: 


s  =  a  exp(i0)  for  1  <  k  <  K;  a  >  0,  0  constant  ..  (67) 


For  these  two  examples,  the  values  of  the  matched  deflection  parameter  in 
(63)  become,  respectively. 


d 


(1) 

m 


^1/2  a 


(68) 


On  the  other  hand,  for  a  clipped  reference,  (65)  yields  the  corresponding 
values 


d 


(1) 

r 


[10.187  a/o  for  K  =  128 
n 


2.613  a/<j 


for  K  =  8 


n 


^1/2  a 


for  all  K  . 


(69) 


Since  d  in  (60)  and  (62)  is  a  complete  descriptor  of  performance  for  the 
linear  input  narrowband  correlator,  it  is  seen  that  for  the  phase -modulated 
tone  of  example  1,  the  clipped  reference  requires 


.(1)\ 


-20  log 


.(1) 

/ 


”-20  log(l0.187AT2F)  =  .91  dB  for  K  =  128~ 
-20  log(2.613/v8’)  =  .69  dB  for  K  =  8 


additional  input  s ignal -to-noi se  ratio  relative  to  the  matched  reference.  On 
the  other  hand,  the  pure  tone  signal  example  2  requires  a  0  dB  difference,  as 
seen  by  reference  to  (68)  and  (69).  In  general,  (63)  and  (65)  reveal  a  dB 
difference,  due  to  clipping  the  reference,  of 


10  log 


I  W' 

_ k _ 

5  -  1  sgr,(s^,)][s^^  t  i 


(70) 


I 


MOMENTS  FOR  NONLINEAR  INPUT  CHANNEL 


In  this  subsection,  the  nonlinearity  g  in  the  input  channel  of  figure  10 
is  present;  that  is,  v^^  is  given  by  (48)  in  terms  of  two  nonlinearities,  one 
each  in  the  in-phase  and  quadrature  channels.  For  a  general  nonlinearity  g, 
the  accumulator  output  is 

^  ■  2  Vk  '  2  “k[5'\r'  '  9(Xk,)]  .  (?') 

k  k 

where  input 

X|^  =  S|^  exp(ie)  =  t|^  nj^  .  (72) 

Signal  ^Sj^j  and  reference  [w^^^  are  also  general  at  this  point. 

Random  variable  e  is  held  fixed  for  the  moment.  Substituting  (72)  in  (71), 


there  follows 
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'  '  5  "kf^'^kr  *■  "kr>  *  ’  «!*ki  *■  "kl']  ■ 


(73) 


We  now  assume  that  K,  the  number  of  samples  accumulated,  is  large,  so 
that  z  is  well  approximated  as  a  complex  Gaussian  random  variable.  In  that 
case,  we  can  concentrate  on  the  two  lowest-order  moments  of  z.  The  p-th 
moment  of  the  real  output  of  one  of  the  nonlinearities  is 


9  (^kr  ^  "kr)  = 


dn 


,,  ,1/2  ,"2\^'"*'kr  * 

(2.)  2= 


=  [ 


dx  0(x)  g^(t^^  4-  <,^x)s  Gp(t^^)  , 


(74) 


where  we  used  the  Gaussian  character  of  the  noise,  (47),  and  (8).  In  a 
simi lar  manner. 


3'(^i  ^"ki^  =%(^ki>  • 


(75) 


The  mean  of  complex  random  variable  z  in  (73)  then  follows  easily  as 


m  =  z 
z 


'5"k[=i(V>  " '  '^i<»ki)]  ■ 


(76) 


To  determine  the  second  central  moments  of  z,  we  combine  (73)  and  (76)  in 
the  form 


z  -  z 


(77) 


where 
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“k  ■  "k&'‘kr  *  "kr>  -  ■ 

“k  ■  '  "k&<‘kl  *  "k1>  -  •  <’ 

From  (47),  every  random  variable  a,  is  independent  of  every  a  for  m  at  k, 

k  m 

and  independent  of  every  for  all  m.  A  similar  property  holds  for  every 

random  variable  B,  .  In  addition.  oT  =  B”  =  0.  It  then  follows  that 
k  k  k 

(z  -7)^.  2  z  »  U  -  I  K  *  8^)- 

km  k  ^  / 


=  I 


W  -  ^('kr'  -  V\i>  *  • 


z 


-2  <» 


1  )  I  +*6 

k  kH  m  m 


)"Z(kI 

k  \ 


•  ^  Kl^[‘52<\r>  -  =?<‘kr>  *  -  '=(<tki)]  ■ 


PROBABILITIES  FOR  CLIPPED-INPUT  CORRELATOR 


We  now  specialize  these  results  to  the  case  of  the  dipper  g(x)  =  sgn(x) 


Then  (74)  yields 


6^(t)  =  j”  5gn(t  o^x)  =  2^(t/o^)  -  1  , 


62(t)  =  J  dx  a(x)  sgn  (t  +  o^x)  =  1  . 


V* 


We  further  restrict  attention  to  the  case  of  small  input  signal -to-noise 


ratio  in  figure  10.  that  is. 


hkl 


«  1  for  all  k  , 


in  which  case  (80)  yields,  with  the  help  of  (32), 


-T-'  S^(tki)  " 


.  II)' hi 


utilization  of  these  approximations  in  (76)  yields  mean 


m  =  - 
z  W 


1/2 


k 


.  t,  +  i  t,  . 
*  kr  ki 


1  /2 

2\  exp( ie) 


k^k 


where  we  used  (72) . 


Similarly,  (79)  reduces  to 


By  the  small  input  signal-to-noise  ratio  assumption  in  (81),  it  can  be 

1  2  _  2 
that  the  magnitude  of  (z  -  z)  in  (84)  is  much  smaller  than  z  -  z 


(81) 


(82) 


(83) 


(84) 

seen 

In 


this  case,  we  have  the  good  approximations 
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(z  -  z)^  s  0  . 

|z  -  Ip  2  2  2. 

k 

Observe,  that  to  this  order  of  approximation,  all  dependence  of  these  moments 
on  signal  has  disappeared. 


Thus,  for  signal  present,  the  joint  probability  density  function  of  the 
real  and  imaginary  parts  of  z  in  (71)  is 


p(Zr.z^) 


2  2 
(y  -  '"zr)  "  (^i  -  ^j) 


(86) 


Then  in  an  identical  manner  to  that  developed  in  (58),  the  detection 
probability  is,  for  threshold  T, 


Pg  =  Prob(-r  >  T)  =  Prob 


=  0(a.  T/oJ  , 


(87) 


>  ".p  -j- r.f  r-f.f  tt, 
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'  j 

.V 

Notice  that  random  variable  e,  which  was  present  in  (72)  and  (83),  has  now 

disappeared  in  A  and  detection  probability  P^.  Also  A  is  precisely  equal  to 

1  /2 

d  in  (60)  except  for  a  scale  factor  of  (2/ir)  .  Thus,  the  clipped-input 

correlator  is  degraded  relative  to  the  linear-input  correlator  by  1.96  dB, 
measured  at  the  system  inputs,  regardless  of  the  signal  [s^^j  and  reference 

Kl- 


If  we  set  the  signal  to  zero  in  (87)  and  (88),  we  obtain  false  alarm 
probabi 1 ity 


P 


F 


0(0,  T/o^)  =  exp/- 


(89) 


where  is  the  common  variable  defined  in  (85)  and  used  in  (86)-(88). 
Eliminating  threshold  T  from  (87)  and  (89),  we  obtain  the  receiver  operating 
characteristic  as 


P 


D 


(-2.?n  Pp)^^^)  . 


(90) 


The  parameter  a  defined  in  (88)  is  a  good  descriptor  of  performance  when  the 
assumptions  of  large  K  and  small  input  signal -to-noise  ratio  are  met. 


Matched  Reference 

The  deflection  parameter  A  in  (88)  is  maximized  by  choosing  a  matched 
reference,  that  is,  w^^  =  ASj^,  thereby  yielding 
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The  last  relation  follows  upon  reference  to  (63).  It  is  seen  that  the 
performance  of  the  clipped-input  narrowband  correlator  depends  only  on  the 
total  received  signal  energy  in  this  case. 

Clipped  Reference 

For  the  clipped  reference,  we  again  make  the  choice  delineated  in  (64). 
The  result  for  the  corresponding  value  of  deflection  parameter  A  in  (88)  is 


where  d^  is  given  by  (65).  Thus,  the  degradation  of  A^,  relative  to 

matched  value  A  above,  depends  on  the  specific  choice  of  signal  waveform 
m 

selected. 


EXAMPLE 


For  the  phase-modulated  tone  presented  in  (66),  (91)  yields 


(93) 


a 

a 


1.96  d8  degradation  relative  to  the  linear-input  narrowband  correlator  with 
matched  reference.  On  the  other  hand,  (92)  and  (69)  yield 

(^8.128  a/o  for  K  =  1281 


A 


(1) 

r 


'2.085  a/o 


n 


for  K 


(94) 


giving  a  loss  that  depends  on  the  particular  number  of  samples,  K. 
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GRAPHICAL  RESULTS 

In  figure  11,  the  receiver  operating  characteristics  (ROC)  for  the 
linear-input  narrowband  correlator,  as  governed  by  (62),  are  drawn  in  dotted 
lines,  for  the  range  (.001,  .8)  in  false  alarm  probability  and  (.001,  .999)  in 
detection  probability.  These  are  no  longer  parallel  straight  lines,  as  they 
were  for  the  baseband  single-channel  case  in  figures  1-9,  but  bend  downward  as 
the  false  alarm  probability  increases  (that  is,  as  the  threshold  decreases). 

Superposed  as  the  jagged  solid  lines  are  the  results  of  two  simulations, 

each  employing  30,000  trials  ,  for  the  linear-input  narrowband  correlator 

with  K  =  128  samples  and  for  the  phase-modulated  tone  example  of  (66).  For 

case  A,  a/o^  =  Y2V4,  whereas  for  case  B,  a/o^  ^  YJ'/B.  These  values  were 

chosen  so  that  the  deflections  would  take  on  the  values  d^^^  ^  4  and 

m 

d^^ ^  =  2,  respectively,  as  may  be  seen  by  reference  to  (68).  The  overlay 
of  simulation  and  theory  is  excellent  in  figure  11  over  the  entire  range  of 
probabilities  plotted,  and  establish  the  degree  of  confidence  to  be  expected 
in  the  simulation  results  to  follow. 

In  figure  12,  K  is  reduced  to  8,  while  the  values  of  a/o^  are  increased 
by  a  factor  of  4,  thereby  realizing  the  same  d-values  as  above,  according  to 
(68).  The  agreement  of  results  is  again  very  good,  except  for  the  upper  end 
of  the  a/o^  -  curve,  where  the  simulation  result  is  low.  This  particular 

*A  sample  program  is  listed  in  the  appendix. 
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simulation  is  atypical;  a  re-run  with  a  different  set  of  30,000  independent 
trials  yielded  good  agreement  over  the  entire  range  plotted.  However,  it  does 
serve  to  illustrate  a  pitfall  of  simulation  results,  even  when  based  on  30,000 
trials,  as  these  results  are;  it  is  possible  to  get  a  systematic  error  in  the 
receiver  operating  characteristic,  especially  on  the  tails  of  the 
distributions. 


In  figure  13,  K  is  increased  back  to  128,  but  the  reference  is  now 
clipped.  According  to  (69),  we  now  have 


,, ,  ,3.60  for  a/o  =  V?/4~) 

d^;)  "  (  for  K  =  128  . 

,1.80  for  a/o„  =  y7/8 


n 


(95) 


J 


which  values  are  borne  out  by  the  simulation  results  in  figure  13. 


When  K  is  decreased  to  8,  (69)  now  yields 


(3.70  for  a/o^  =  )[? 

^  ll  .85  for  a/o  =  ^/{T 
^  n 


for  K  =  8 


(96) 


These  results  are  confirmed  by  the  plots  in  figure  14. 


When  the  input  to  the  narrowband  correlator  is  clipped,  instead  of  the 
reference  (which  is  matched  according  to  (91)),  the  pertinent  equation  is 
(93) .  There  fol lows 

(*3. 19  for  a/o  ■  (7/4^ 

a;  ^  f  "  (  for  K  -  l?S  .  (9/) 

( 1 .60  for  a/o  ^7/ti 
n 


figure  15  corroborates  these  predictions. 
46 
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1 

When  K  is  decreased  to  8,  while  a/o  is  increased  by  a  factor  of  4,  the 
same  values  result  from  the  use  of  (93),  as  given  in  (97).  The  simulation 

results  in  figure  16  reveal  that  the  actual  performance  of  the  narrowband 
correlator  does  not  meet  these  deflection  values  3.19  and  1.60,  except  at  the 
upper  right  end  of  the  curves.  The  reasons  for  the  discrepancy  are  two-fold: 

K  =  8  is  not  large  enough  to  rely  on  the  Gaussian  approximation,  and  input 
signal-to-noise  ratios,  a/o^  =  'V?  and  l/vT,  are  not  small,  as  assumed  in 
(81 )  et  seq. 


When  both  the  input  and  the  reference  are  clipped,  the  simulation  results 
for  K  =  128  are  displayed  in  figure  17.  The  corresponding  theory  is  furnished 
by  (92)  and  (65),  and  yields  the  results  for  already  listed  in 

(94) .  They  become 

(2.87  for  a/0  = 

i  "  >  for  K  =  1  28  ,  (98) 

(1.44  for  a/o^  =  V?/8  J 

and  are  1n  excellent  agreement  with  the  simulations  In  figure  17. 


For  the  alternative  situation  with  K  =  8,  the  appropriate  values  follow 
from  (94)  as 


,(1)  ,[^  ’5  'O'- 

1  .47  for  a/o  =  1/V2  > 
n 


for  K  ^  8  . 


(99) 


These  values  furnish  a  good  approximation  to  the  simulation  results  in 
figure  18  at  the  upper  right  end,  but  are  optimistic  over  the  rest  of  the 
range.  The  reasons  for  the  discrepancy  are  again  that  K  is  not  large  and  the 
input  signal-to-noise  ratio,  a/o  ,  is  not  small. 


001  L-J  ■] _ ^ ^ _ ^ _ L__J _ ^ _ ] _ 

.001  .01  .02  .05  .1  .2  .3  .4  .5  .G  .7 

Probability  of  Faloo  HI  arm 


Figure  16.  ROC  for  Narrowband,  Clip  Input 
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Figure  18.  ROC  for  Nar  roujb  and  ,  Clip  Both 
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The  final  result  In  figure  19  Is  Identical  to  the  conditions  of  figure 

17,  except  that  now,  the  quadrature  component  of  the  clipped  reference 

channels  was  suppressed;  that  is,  only  one  real  reference  was  used  for 

multiplication  in  figure  10.  The  simulation  results  in  figure  19  indicate 

values  of  approximately  2.6  and  1.3;  these  values  are  approximately 
r 

.9  dB  poorer  than  (98)  which  pertained  to  figure  17.  Thus,  dropping  one  of 
the  clipped  reference  channels  causes  a  loss  of  almost  1  dB,  and  should  be 


avoided . 


gs 


DISCUSSION/SUMMARY 


Accurate  prediction  of  performance  of  the  baseband  and  narrowband 
correlators,  with  clipping  at  the  input  and/or  reference  levels,  is  possible 
over  a  wide  range  of  number  of  input  samples,  input  signal-to-noise  ratio,  and 
detection  and  false  alarm  probabilities.  Only  when  K,  the  number  of  samples, 
gets  too  small  to  use  the  central  limit  theorem,  does  the  theory  begin  to 
deviate  from  actual  performance.  However,  it  is  precisely  in  this  case  that 
simulation  is  most  attractive,  since  a  large  number  of  trials  can  then  be 
conducted  in  a  reasonable  amount  of  time.  Programs  are  furnished  in  the 
appendix,  for  both  the  baseband  as  well  as  the  narrowband  correlator,  that 
enable  extension  of  the  simulation  results  to  other  cases  of  interest  to  the 
user. 


If  we  compare  corresponding  results  in  the  baseband  and  narrowband 
correlators,  that  is,  figure  2  with  figure  11,  figure  3  with  figure  12,  etc., 
the  predicted  and  simulation  results  for  the  deflections  are  identical,  with 
one  exception.  Namely,  when  the  reference  is  clipped  and  K  is  small,  the 
narrowband  deflections  are  slightly  larger  than  the  corresponding  values  for 
the  baseband  correlator;  compare  figure  5  with  figure  14,  and  compare  figure  9 
with  figure  18.  Thus,  the  degradation  suffered  in  the  case  of  a  narrowband 
correlator  is  not  quite  as  bad  as  for  the  baseband  correlator  in  these 
particular  cases. 
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AfffNDIX 


PROGRAMS 

The  first  program  listed  in  this  appendix  pertains  to  the  baseband 
correlator  of  figure  1,  with  both  the  input  and  reference  clipped;  see  lines 
320-350.  If  clipping  is  desired  removed,  merely  delete  SGN  from  the 
appropriate  lines.  This  program  is  written  in  BASIC  for  the  Hewlett-Packard 
9000  computer;  the  qualifier  DOUBLE  in  line  60  denotes  INTEGER  variables. 

The  second  program  pertains  to  the  narrowband  correlator  depicted  in 
figure  10,  with  both  the  input  and  reference  clipped;  see  lines  370-400.  If 
clipping  is  desired  removed,  delete  the  appropriate  SGN  operations. 
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1  0 

'  CLIP  INPUT  FIND 

REFERENC 

2  0 

I  =  3  0  0  0  0 

1 

3  0 

I  =8 

1 

4  0 

Fl  =  =2 . 

1 

^■0 

RT="BEE2" 

60 

DOUBLE  I,K,Ks. , 

I  5-  ' 

70 

D  I  N  Z  '■  I  ;  80000  > 

, Cos  *1:1 2 

S  0 

REDIN  Z‘.1:I  ',C 

‘  1 : K > 

■30 

P'4  =  -L0G':4.  > 

1  0  0 

T  =  2.  *PI  .  K 

1  10 

FOR  Ki=l  TO  y 

no 

F‘  h  1  k  =  T  *  K  = 

1  30 
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1  ' 
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NEXT  Ki 

1  ^'0 
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j 
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1 0  0 
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I  CLIP  INPUT  hND  PEFEPENCE; 

1  =  J  0  0  01  01  I 

C  =S  ' 

P)  £  =  p  C'  P  ‘  2  .  '  ' 
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